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Abstract: The mitogen-activated protein kinase-activated protein kinase MK5 is a 
substrate of the mitogen-activated protein kinases p38, ERK3 and ERK4. Cell culture and 
animal studies have demonstrated that MK5 is involved in tumour suppression and 
promotion, embryogenesis, anxiety, cell motility and cell cycle regulation. In the present 
study, homology models of MK5 were used for molecular dynamics (MD) simulations of: 
(1) MK5 alone; (2) MK5 in complex with an inhibitor; and (3) MK5 in complex with the 
interaction partner p38a. The calculations showed that the inhibitor occupied the active site 
and disrupted the intramolecular network of amino acids. However, intramolecular 
interactions consistent with an inactive protein kinase fold were not formed. MD with p38a 
showed that not only the p38 docking region, but also amino acids in the activation 
segment, aH helix, P-loop, regulatory phosphorylation region and the C-terminal of MK5 
may be involved in forming a very stable MK5-p38a complex, and that p38a binding 
decreases the residual fluctuation of the MK5 model. Electrostatic Potential Surface (EPS) 
calculations of MK5 and p38a showed that electrostatic interactions are important for 
recognition and binding. 
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1. Introduction 

Mammalian mitogen-activated protein kinases (MAPKs) are cytoplasmic protein- Ser/Thr kinases 
that participate in signal transduction by catalysing the transfer of the y-phosphoryl group from ATP to 
a hydroxyl group of the protein substrate. All mammalian cells possess multiple MAPKs involved in 
controlling cellular processes such as proliferation, differentiation, survival, apoptosis, gene regulation 
and motility [1]. 

To date, seven distinct MAPKs pathways have been discovered in mammalian cells, and the 
MAPKs are divided into conventional and atypical MAPKs. The conventional MAPKs constitute three 
consecutive phosphorylation events mediated by three Ser/Thr protein kinases. This three-tiered 
cascade consists of: (1) MAPKs; (2) MAPK kinases that phosphorylate and activate MAPKs; and 
(3) MAPK kinases kinases that phosphorylate and activate MAPK kinases. Conventional MAPKs 
include the extracellular signal-regulated kinases 1 and 2 (ERKl/2), p38 MAPKs (p38 a, (3, y, and 5), 
c-Jun terminal kinases 1-3 (JNKl-3) and the ERK5 [1]. Conventional MAPKs are characterized by a 
conserved Tyr-X-Tyr motif in the activation loop of the kinase domain. Atypical MAPKs lack this 
motif, and are not organized in the classical consecutive phosphorylation cascade of three Ser/Thr 
protein kinases. Atypical MAPKs comprise ERK3/4, ERK7/8 and the Nemo-like kinases (NLK) [2]. 

Activated MAPKs catalyse the phosphorylation of a large number of proteins including protein 
kinases, phosphatases, transcription factors and several other proteins. Protein kinases that are 
phosphorylated by the MAPKs are termed MAPK-activated protein kinases (MAPKAPKs), and so far 
1 1 human Ser/Thr MAPKAPKs have been identified. Based on sequence similarities, MAPKAPKs are 
subdivided into four groups: (1) Ribosomal S6 kinases (RSKl-4); (2) Mitogen and stress-activated 
kinases (MSKl and 2); (3) MAPK-interacting kinases (MNKl and 2); and (4) MAPKAPK 2, 3 and 5 
(MK2, MK3, and MK5) [3,4]. 

All MAPKAPKs are phosphorylated and activated by different conventional MAPKs. However, 
MK5 is different from the others since MK5 is phosphorylated and activated both by conventional 
(p38) and atypical (ERK3 and ERK4) MAPKs (Figure 1) [1,2,5]. ERK3 and ERK4 have been 
suggested to catalyse phosphorylation at Thrl82 and activate MK5 [6,7]. 

It is generally accepted that MK2 and MK3 are bona fide p38 substrates [4,8]. However, the in vivo 
interaction between p38 and MK5 is under some debate, and is currently not completely resolved 
(reviewed in [9]). Already in 1998 it was demonstrated that recombinant p38 could phosphorylate 
MK5 in vitro [10]. Thrl82 of recombinant MK5 was also found to be phosphorylated by p38a, p38p, 
p385 and p38y in vitro, but only p38a and p38p enhanced the intrinsic enzymatic activity of MK5 
towards the recombinant small heat shock protein HSP27 [11]. Furthermore, it was demonstrated that 
p38 regulates the subcellular localization of MK5. Endogenous MK5 and EGFP- or HA-tagged MK5 
are predominantly located in the nucleus, but after stress stimulation by TNFa, sorbitol, and sodium 
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arsenite, the p38 pathway is activated and causes nuclear export of MK5 [12,13]. Moreover, 
experimental studies indicated that nucleocytoplasmic shuttling is controlled through interaction with 
PKA, CdclSA and ERK3, ERK4 and p38 [7,12,14-19]. However, yeast two-hybrid assays showed 
that MK5 interacts only weakly with p38a, and not at all with p385 and p38y [6]. Another yeast 
two-hybrid study identified MK2 as a p38 interaction partner, but not MK5 [7]. 



Figure 1. MAPK signalling cascades leading to MK5 activation. Mitogen-activated protein 
kinase kinase kinases 1-4 (MEKKl-4), the dual leucine zipper kinase (DLK), the mixed 
lineage kinase (MLK2), the mitogen-activated protein kinase kinase kinase 8 (Tpl-2), the 
apoptosis signal-regulating kinase 1 (ASKl), the transforming growth factor-(3-activated 
protein kinase (TAKl) and the serine/threonine-protein kinase (TAOl/2) phosphorylates 
and activates dual specificity mitogen-activated protein kinase kinase3/6 (MEK3/6), which 
in turn phosphorylates and activates p38 a, (3, 5 and y. p38 will then phosphorylate and 
activate mitogen- and stress-activated protein kinases 1/2 (MSKl/2), the MAPK-interacting 
kinase 1 (MNKl), MAPK-activated protein kinases 2/3 (MK2/3), and MAPK-activated 
protein kinase 5 (MK5). In addition to the p38 pathway MK5 can also be activated by 
extracellular signal-regulated kinases 3/4 (ERK3/4). ERK3/4 are activated by p21 -activated 
kinases 1/2/3 (PAK 1/2/3). The kinase/kinases activating PAK is currently unknown. 
Dotted line indicates that although activation has been reported, it has to be thoroughly 
investigated, while solid lines indicate validated activations. 
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MK5 is expressed in all examined cell types and tissues, but most abundantly in heart, brain and 
hematopoietic progenitors [11,20-22]. Several experimental studies have been performed to elucidate 
the biological role of MK5. These studies have shown that MK5 may be involved in a wide range of 
biological processes including control of the cytoskeletal architecture by F-actin remodelling [12,23-26], 
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tumour growth and progression by stimulating angiogenesis [27], tumour suppression by Myc 
downregulation or antagonizing the JNK pathway [28,29], or interfering with ERK3-mediated cell 
cycle control, and B-cell development by activating Rag transcription via Foxol [18,30-32]. Further, 
MK5 was also suggested to be involved in the inflammatory response as a consequence of microbial 
infection [33] and heart function [2]. Animal studies have also suggested a role of MK5 in neurological 
processes controlling anxiety and locomotion [34,35]. However, the exact biological role of MK5 is 
not completely understood, and the main reasons are lack of selective MK5 inhibitors and that MK5 
knock-out mice do not display any obvious phenotype under normal conditions [7,36,37]. 

Structural information about the MK5 and its molecular interactions that can be used to design 
selective MK5 compounds is therefore of great interest. Regulation of MK5 activity by compounds 
may also be therapeutically valuable in different disease conditions. The exact three-dimensional 
structure of MK5 has not been resolved by X-ray crystallography. However, the X-ray crystal 
structures of MK2 in complex with p38a [38], ADP [39] or small molecular inhibitors [40-43], and of 
MK3 in complex with an inhibitor [44] are known. Because MK5 has an overall sequence identity of 
42% with MK2, and 41% with MK3, it is believed that MK5 may have a 3D fold similar to that of 
MK2and MK3. 

We have previously constructed 3D homology models of the human analogue of MK5. Our 
calculations indicated that the models were of high quality and had predictive values [45]. The MK5 
models display the protein kinase fold which is well conserved among all members of the protein 
kinase family. In the present paper, we studied the molecular interactions of MK5 with interaction 
partners. The program Desmond, as implemented in the Schrodinger suite of programs, was used for 
comparative molecular dynamics (MD) simulations of MK5. The previously published MK5 
homology models were used as starting structures for MD simulations with MK5 alone, in complex 
with a pyrrolopyrimidone-based inhibitor and in complex with p38a. 

2. Results and Discussion 

2. 7. The Protein Kinase Fold and Homology Models of MK5 

The typical bilobal protein kinase fold consists of a small A^-terminal lobe and large C-terminal 
lobe connected by a structurally flexible linker (hinge) region. Our MK5 models show that this fold 
is also conserved in MK5 (Figure 2). In MK5, the hinge region corresponds to the segment 
Glul03-Metl04-Metl05. The residue A^- terminal of the hinge region is termed the gatekeeper residue 
(Met 102 in MK5). The X-ray structure of ERK2 shows that the gatekeeper residue confers selectivity 
for binding nucleotides and small molecular inhibitors, but also seems to control auto-phosphorylation 
through intramolecular interactions [46]. The A^- terminal lobe of ERKl and ERK2 includes a 
five-stranded antiparallel P-sheet (pi-p5) and a highly conserved aC-helix. The C-terminal lobe is 
mainly a-helical and consists of six conserved a-helixes (aD-al) and four small P-strands [47]. 
Figure 2 is showing that the N- and C-terminal lobes of the MK5 model have an overall organization 
similar to that of ERKl and ERK2. 

X-ray structures of protein kinases show that the active site is located in between the two lobes and 
constitutes the binding of ATP and two magnesium ions. The active site is shielded by a loop named 
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the ATP-phosphate binding loop (also called the P-loop). The P-loop is located between (31 and (32 in 
the N-lobe (Figure 2). The loop contains a glycine-rich segment (GxGxxG), and in ERKl and ERK2 
this loop is suggested to be important for positioning the ATP (3- and y-phosphates for phosphate 
transfer to the substrate [47]. The GxGxxG domain is also conserved in MK5 and corresponds to the 
Gly29-Ala30-Gly31-Ile32-Ser33-Gly34. The GAGISG domain is also conserved in MK5 sequences 
from other species [25]. 

Figure 2. The Ca-trace of the MK5 homology model based on the X-ray structure of MK3 
in complex with an inhibitor (PDB id: 3FHR). The construction of the model was 
described in Lindin et al. [45]. 
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X-ray structures of ERKl and ERK2 have shown that the aC-helix of the N-lobe may occur in an 
activated and inactivated orientation. Structures of inactive and active ERKl and ERK2 have shown 
that in the activated structure a glutamic acid of aC forms a salt bridge with the lysine of an AXK 
sequence in the (33 strand. In MK5, the aC helix consists of residues Pro57-Ala71, and the glutamic acid 
corresponds to Glu62. The (33 strand is constituted by Arg47-Leu54 in MK5 and the lysine corresponds to 
the Lys51 of the Ala49-Leu50-Lys51 sequence. The active conformation is often named the aC-in 
conformation, while the inactive conformation is named the aC-out conformation of protein kinases. 

Three amino acids in ERKl and ERK2 are defining a catalytically important K/D/D motif [47]. 
These amino acids are also conserved in MK5. In ERKl and ERK2, the lysine of this motif resides in 
the (33 strand and corresponds to Lys51 in MK5. In addition to forming a salt bridge with the glutamic 
acid in aC in the activated state, the lysine also binds the a- and (3-phosphorus atoms of ATP. The 
aspartic acids of the K/D/D motif are located in the C-terminal lobe as part of the catalytic loop and the 
activation segment, respectively (Figure 2). In MK5, the first aspartic acid of this conserved K/D/D 
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motif is part of a HRDLK segment in the catalytic loop, located in the C-terminal of the p6-strand, and 
corresponding to Hisl46-Argl47-Aspl48-Leul49-Lysl50-Prol51 in MK5. In both ERKl and ERK2 
this aspartic acid is suggested to position the substrate hydroxyl group, and extract the proton of the 
OH-group and facilitate a nucleophile attack of the substrate oxygen on the y-phosphorus atom of 
ATP. The aspartic residue thereby functions as a base in the catalytic reaction [47]. The ERKl and 
ERK2 structures also show that the lysine (corresponding to LyslSO in MK5) of this segment binds the 
oxygen of the substrate hydroxyl group and the y-phosphorous atom of ATP. The second aspartic acid 
in the K/D/D motif of ERKl and ERK2 is the first amino acid in the activation segment that begins 
with the DFG motif corresponding to Aspl69-Phel70-Glyl71 in MK5. The DFG motif is followed by 
the activation loop that contains the regulatory phosphorylation site (Thrl82 in MK5). In an active 
protein kinase, the aspartic acid is pointing into the active site and interacts with two Mg ions that 
interact with the a-, P- and 6 phosphates of ATP. The activation loop is fundamental both for substrate 
binding and catalysis. 

In our previous study [45] we constructed several homology models of MK5. In the present study 
these models were used for MD simulations. The model constructed by using the X-ray structure of 
MK3 in complex with an inhibitor [44] as template, showed highest accuracy in selecting known MK5 
binders in front of decoys during docking. The model based on the MK3 structure was therefore used 
for 100 ns MD simulation of MK5 alone and for two simulations (100 and 200 ns) of MK5 in complex 
with a pyrrolopyrimidone-based inhibitor. The MK5 model used for MD with p38a was based on the 
X-ray crystal structure of MK2 in complex with p38a [38,45]. 

2.2. Stability of the Molecular Systems 

High quality X-ray structures are first choice as starting structures for MD simulations and other 
predictions of molecular mechanisms of action using molecular modelling. However, in the absence of 
such structures, simulation with high quality homology models may contribute with important 
structural and mechanistic insight into different classes of proteins and their interactions. Proteins and 
protein mechanisms that recently have been successfully studied by combining homology modelling 
and MD simulations include among others: the E7 protein from high- and low-risk types of human 
papillomavirus [48], the CHeW coupling protein of the chemo taxis signalling complex [49], lipases 
and their interactions with known inhibitors [50], molecular mechanism of the serotonin transporter [51], 
conformational changes underlying ion channel opening [52] and molecular mechanism of 
G-protein-coupled receptors [53-55]. 

After using the default relaxation/equilibration protocol of the Desmond program, which includes 
series of steepest decent energy minimisations and short MD simulations at different temperatures 
(described in Section 3.3), the molecular systems were simulated for 100 ns. The plots of Root Mean 
Square Deviations (RMSD) (Figure 3), potential energy (Figure 4) and Root Mean Square Fluctuation 
(RMSF) (Figure 5) and atomic distance during 100 ns sampling simulations do not include the 
relaxation/equilibration period. The RMSD of backbone atoms of the MK5 models relative to initial 
structures show that stable molecular systems were obtained during MD. The RMSD plots indicated 
that the free MK5 model and the MK5-inhbitor complex were more unstable in the first part of the 
simulations than the MK5-p38a complex. The template for constructing the MK5 model that was 
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used for MD of the free enzyme and the MK5 -inhibitor complexes was a MK3 -inhibitor complex 
(PDB id: 3FHR), and some structural reorganization may be necessary to account for the removal of 
the MK3 inhibitor during the modelling process. During 100 ns MD with the p38a complex the RMSD 
of the backbone stabilised at around 3 A after approximately 30 ns of simulation. During 100 ns MD 
with the inhibitor, the complex stabilised at around 7.5 A after about 45 ns, while the free MK5 model 
stabilised at around 5 A after about 45 ns of MD (Figure 3). The trend for the 200 ns MD with the 
inhibitor was similar to the 100 ns MD (data not shown). 

Figure 3. The Root Mean Square Deviations (RMSD) of backbone atoms relative to the 
starting complexes during 100 ns MD. 
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Figure 4. Potential energy of the molecular systems during 100 ns MD simulation. 
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The potential energy of the system as the sum of bond, angle, torsion, and non-bonded terms can 
also be used as a simple measure of system stability. Thus, plots of potential energy as a function of 
time were generated for the simulated complexes (Figure 4). The plots indicated that all molecular 
systems were well equilibrated and remained stable throughout 100 ns simulations. This was also 
observed during a 200 ns MD with the inhibitor (data not shown). The RMSD and the potential 
energies show that we have obtained structurally and energetically robust models. 

Figure 5. The Root Mean Square Fluctuation (RMSF) during 100 ns MD. All atoms were 
included in the RMSF calculations. 
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2.3. MD of the Free MK5 Model 

The structurally most flexible regions during 100 ns MD of the native MK5 model were as expected 
the loop regions, and particularly the region around residues 190-210 (Figure 5). This region is at the 
C-terminal end of the activation loop. In most protein kinases the activation segment starts with the 
DFG motif (Asp 169-Phel70-Gly 171 in MK5) and ends with an APE motif. The APE motif is present 
in both MK2 and MK3, but in MK5 this motif has been replaced with an APQ motif. The APQ motif 
seems to be conserved between MK5 from different species. Further, the MK3 structure (PDB: 3FHR) 
that was the template for constructing the MK5 model is disordered in this part of the activation loop, 
which may indicate that this part of the activation loop is also structurally flexible in MK3. 

In addition, the MK5 has an insert between the APQ motif and the aF-helix compared with MK2 
resulting in a longer activation loop than in MK2 and MK3. The lacking part of this loop was 
constructed by homology with available loop structures in the PDB database [45], and is present at the 
surface of the MK5 model. 

The same MK5 model was used as starting structure both for the MD of MK5 alone and for the 
MD of the MK5-inhbitor complex (Figure 6). The starting MK5 model was in an active kinase 
conformation with interactions between Lys51 and Glu62 (Figures 6 and 7), while Asp 169 of the DFG 
motif in the activation segment interacts with amino acids in the P-loop. The interacting distance 
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between the side chain terminal atoms of Lys51 and Glu62 was varying between 2.5 or 4.0 A up to 
10 ns, then the interacting distance was in the range of 2.5-6 A, and that lasted up to about 30 ns 
(Figure 7). After 30 ns, the distance stabilised and was flipping between the range of 1.5-2.3 A and the 
range of 2.8-4.0 A and that lasted up to the end of the simulation. Further, Lys51 also interacted with 
Ser33 (Figure 6) located in the P-loop. In the beginning of the simulation, the binding distance 
between the hydroxyl group of Ser33 and Lys51 was around 4 A, but at the same time as the 
Lys51-Glu62 atomic distance stabilised (at 30 ns), the Lys51-Ser33 atomic distance (Figure 6) 
stabilised and was flipping between the range of 1.5-2.3 A and the range of 2.8-4.0 A, similar to the 
Lys51-Glu62 atomic distance. After 30 ns, we observed concerted changes in the interacting distances 
of the amino acid triplet Ser33, Lys51 and Glu62 (Figure 6). When the atomic distance of Ser33-Lys51 
was in the shortest range (1.5-2.3 A), the atomic distance between Lys51 and Glu62 was also in the 
shortest range (1.5-2.3 A), and when Ser33-Lys51 was in the longest range, the Lys51-Glu62 distance 
was also in the longest range. Glu62 also interacted with backbone NH-group of Phel70 (DFG motif). 
In the beginning of the MD, the atomic distance varied between 2.0-4 A. However, at 30 ns when the 
interactions of Glu62-Lys51 and Lys51-Ser33 stabilised, the interactions between Glu62 and PhelTO 
also stabilised, and stayed around 2.0 A for the rest of the MD. Asp 169 in the DFG motif, which is 

2_|_ 

coordinating the Mg ions and the a and y phosphate groups when ATP binds, interacted with Ile32 
and Ser33 during MD (Figure 6). Asp 169 interacted both with the backbone NH-group and the side 
chain hydroxyl group of Ser33. The Aspl69-Ser33 backbone interactions varied in the range of 
2.5-4.5 A up to 33 ns, and then stabilised in the range 1.7-2.3 A, while Aspl69-Ser33 side chain 
interactions varied between 3-4.5 A up to 33 ns and then both stabilised in the range of 1.7-2.3 A for 
the rest of the MD. The interactions with the backbone NH-group of Ile32 were in the range of 2-3.5 A 
up to about 43 ns, but then stabilised at about 2 A. 

Figure 6. The network of interacting amino acids at the active site of the starting structure 
for MD with MK5 alone and in complex with the inhibitor. The model was based on an 
X-ray structure of MK3 in an active kinase conformation, and the MK5 model is therefore 
also in active kinase conformation. Colour coding of atoms: red; oxygen, blue; nitrogen, 
green; carbon atoms of the inhibitor, grey; carbon atoms of the amino acids, white; 
hydrogen atoms. Dotted lines are showing observed interactions between amino acids 
during MD of the free MK5 and the MK5 -inhibitor complex. 
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Figure 7. The atomic distances (A) between Lys51 and Glu62 during 100 ns MD 
simulations, (a) free MK5; (b) inhibitor-MK5 complex; (c) p38a-MK5 complex. 
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The intermolecular network of interactions, including Ser33, Lys51, Glu62 and Phi 70, stabilised at 
around 30 ns, while the interactions of Asp 169 stabilised from between 30 and 40 ns of MD. For the 
rest of the MD, the interactions of Aspl69 with Ser33 and the interactions of Phel70 with Glu62 were 
in the region of 1.7-2.5 A, while the Lys51-Ser33 and Lys51-Glu62 interactions were flipping between 
1.5-2.3 and 2.8-4.0 A. Taken together, the MD simulation of the free MK5 model indicated that 
the hydrogen bonding network between amino acids at the active site, consistent with an active open 
protein kinase fold, was stably formed after about 30 ns of MD and was preserved throughout the MD. 

2.4. The MK5 -Inhibitor Complex 

Incorrect ligand protonation may influence ligand-protein interactions during MD. However, 
predicting the pKa value is not straightforward, especially when several ionisable groups exist in 
the ligand. According to CHEMBEL database (https://www.ebi.ac.uk/chembl/), the pKa of the 
pyrrolopyrimidone-based inhibitor is 6.58. pKa predictions using the Jaguar program of the 
Schrodinger suite of programs, and the ICM-software indicated a pKa value in the range of 6.5-6.6. It 
is therefore highly possible that the majority of the inhibitor molecules are in neutral form when 
interacting with MK5. However, it is still quite possible that some MK5 molecules may interact in 
protonated form. The most probable group for being protonated is the nitrogen atom of morpholine 
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ring. Docking of both the neutral (Figures 6 and 8) and protonated forms of the inhibitor showed that 
the morphoUne ring did not participate in direct interactions with MK5, but instead interacts with water 
molecules on the surface of MK5 (Figure 8). Based on this, we decided to do MD simulations with a 
neutral form of the inhibitor. 

Figure 8. 2D plot of the stable inhibitor binding pose observed during the 100 ns MD. 
Amino acids and water molecules within 4 A from the inhibitor are included in the figure. 
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The MD simulation showed that the loops were also the most flexible regions during 100 ns MD 
with the inhibitor, and the region around residues 190-210 showed the largest fluctuation (Figure 5). 
Another very flexible segment was Pro90-Arg95, which was even more flexible than during MD with 
the free MK5 model. This is a part of the loop between p4 and (35 (Figure 2). In MK5, this loop has an 
insert of 6 amino acids compared to MK2 and MK3. Further, the region around residue 240-260 was 
also more flexible during MD with the inhibitor than during MDs with the free model, and with p38a. 
This loop structure is a part of the activation segment and has insertions in MK5 compared to MK2 
and MK3. 

The pyrrollopyrimidone-based inhibitor (Figure 8) was first discovered as a MK2 inhibitor, but was 
later shown to inhibit MK5 and other protein kinases [56]. The MD simulation indicated that binding 
of the inhibitor strongly affected the hydrogen bonding network and the dynamics of the active site. 
The interactions between Asp 169 (DFG motif) and ne32 and Ser33 of the P-loop (Figure 6), that were 
observed for the free enzyme, were not present during 100 ns MD with inhibitor (Figure 9). The 
atomic distance between Glu62 and the backbone NH-group of Phel70 that was around 2 A for most 
of the MD of the free enzyme was also disturbed by the binding of the inhibitor. The Glu62-Phel70 
atomic distance (Figure 6) was around 4 A up to 40 ns, then increased to 6 A for the period between 
40 and 65 ns, and after 65 ns the atomic distance was around 4 A (Figure 9). During the MD with 
the free enzyme, Lys51 formed stable interactions both with Glu62 and with Ser33 in the P-loop 
(Figure 6). The Lys51-Glu62 interactions were also observed for most of the 100 ns simulation with 
the inhibitor (Figure 7); however, close binding interactions between Lys51 and Ser33 were not 
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observed (Figure 9). During the 100 ns of MD with inhibitor, a structure with Lys51 in (33 and Glu62 in 
aC in close proximity, similar to an active protein kinase fold, was obtained, and the interacting 
distance was varying between 3 and 4 A for the entire simulation (Figure 7). 

Figure 9. Atom distances between central amino acids at the binding site during 100 ns 
MD with the inhibitor. 
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The 100 ns MD with the inhibitor indicated that a very stable MK5-inhbitor complex was formed 
and that the inhibitor was in a stable binding mode throughout the MD. During the 100 ns, the 
pyrimidine ring system interacted within a charged and polar part of the MK5 binding pocket, and the 
CO-group of the pyrimidine ring interacted with Lys51 in the aC helix (Figure 8). Up to 25 ns, the 
interacting distance was in the range of 4-6 A, but from 25 ns, stable hydrogen bonding interactions 
(1.8-2.3 A) were observed. The nitrogen atom of the pyrimidine ring interacted with Asnl53, and from 
25 ns to the end of the simulation the interacting distance was fluctuating between 1.7 to 4 A. The 
pyridine ring interacted in a hydrophobic region of MK5 and the nitrogen interacted with Metl05 in 
the hinge region. Up to 35 ns, the interacting distance was about 2 A, and after 35 ns up to the end of 
the MD the interacting distance was fluctuating between 2 and 5 A. The morpholine ring and the 
phenyl ring interacted with amino acids in the C-terminal end of (31 (Gln26-Lys27-Leu28). However, 
water molecules also seem very important for binding the inhibitor (Figure 8). During MD, water 
molecules were forming hydrogen bonds with both the oxygen and the nitrogen atom of the 
morpholine ring, the NH-group of the pyrrolo ring and the NH-group of the pyrimidine ring (Figure 8). 

In order to further test the stability of the inhibitor binding mode, we started a new MD simulation 
of the MK5 -inhibitor complex, but now for 200 ns. The 200 ns MD was started from the same 
equilibrated complex as the 100 ns MD. During the 200 ns, the inhibitor binding mode (data not 
shown) was similar to the binding mode during the 100 ns. The atomic distances between central 
amino acids at the active site were also similar to that observed during the 100 ns MD. 

Taken together, the simulations with the inhibitor indicated that the inhibitor binds to an active 
protein kinase fold of MK5 and functions as an inhibitor by occupying the ATP binding site. Upon 
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inhibitor binding, the intramolecular network and dynamical motions of amino acids at the active site 
are affected (Figures 6 and 8). The intramolecular interactions involving Phel70 (DFG motif) were 
changed, but the interactions between Lys51 and Glu62 were maintained, indicating that a fully 
inactive kinase fold was not formed during the MD. 

2.5. TheMK5-p38 Complex 

The MK5 protein includes functional binding motifs for the kinases ERK3, ERK4 and p38. Further, 
experimental studies have indicated that p38 can phosphorylate and activate MK5 and regulate the 
subcellular location of MK5 [10,11,14,15]. The interactions have been demonstrated in vitro after 
overexpression of MK5 and p38, but the in vivo value of these observations is still under some debate, 
as previously explained, and is not completely resolved [9]. In resting cells, MK5 is predominantly 
located in the nucleus but is able to shuttle between the nucleus and the cytoplasm. MK5 contains 
functional nuclear export signals (NES) and nuclear localization signals (NLS), and the opposite action 
of these motifs may explain the nucleocytoplasmic shuttling of MK5. Activation of the p38 pathway 
was shown to induce nuclear export of MK5. Both p38a and p38p were reported to control the distinct 
subcellular localization of MK5 [14,15,57]. Another goal of the present study was therefore to study 
the molecular interactions between MK5 and p38a. 

The MK2 structure in the MK2-p38a X-ray complex (PDB id: 20ZA) was used as a template for 
constructing the MK5 model of the MK5-p38a complex. The MK5 sequence consists of 473 amino 
acid residues, while MK2 consists of 400 amino acids, such that based on the template, structure we 
could not predict the structure of the entire MK5. The MK5 model is therefore C-terminally truncated 
and consists of the amino acids Met7-Gly367. Previous experimental studies have shown that the 
region Asn356 to Ser373 of MK5 is important for interacting with p38 and was termed the p38 
docking site [14]. This indicates that the model did not contain the entire experimentally verified p38 
docking site. However, experimental studies have also shown that a C-terminal truncated form of MK5 
consisting of residues 1 to 368 also binds p38a, although not as strong as the full length MK5 [14]. 
The docking of p38a into MK5 was guided by the X-ray structure complex of MK2 and p38a. 
The MK5-p38a complex is depicted in Figure 10. The docked complex showed that not only 
Asn356-Gly367 in the C-terminal of the model is involved in p38a binding, but also other amino acids 
in the C-terminal region are also required (Table 1). In addition, the complex showed that amino acids 
in the P-loop, activation segment, aH helix, and the regulatory phosphorylation region are also 
important for binding to p38a (Table 1). The complex showed that MK5 and p38a interact such that 
the MK5 N-lobe interacts with the p38a N-lobe, and the MK5 C-lobe interacts with p38a C-lobe 
(Figure 10). The docked complex indicates that the C-terminal segment which contains most of the 
experimentally verified p38a docking site and the regulatory phosphorylation region of MK5 interact 
(as two arms) at opposite sides of the p38a C-lobe (Figures 10 and 1 IC). 
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Figure 10. The model of MK5 (green) in complex with p38 (yellow), outlining the 
important structural elements of MK5 interacting with p38. The Ca-traces of the proteins 
are shown. 
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Figure 11. Skin representation of the molecular electrostatic potential surface (EPS) 
coloured according to the electrostatic potentials. Colour coding: red; ESP < -5 kcal/e.u. 
charge units, blue; >+5 kcal/e.u. charge units, grey: between -5 and +5 kcal/e.u. charge 
units. (A) Entire complex in the same view as in Figure 8; (B) Entire complex rotated 
180 degrees (y-axis) compared with (A); (C) The EPS of MK5 with backbone Ca-trace of 
p38a in yellow. Same view as in (A); (D) Same as in c, but rotated 90 degrees (y-axis) 
compared with (C); (E) The EPS of p38a and the backbone Ca-trace of MK5 in green. The 
same view as in (A); (F) Same as in (E), but rotated 90 degrees (y-axis) compared 
with (E). 
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Figure 11. Cont. 
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Table 1. Interacting pair of amino acids forming at the MK5-p38a interface during 100 ns MD. 
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Figure 12. A close up of interacting amino acids forming the MK5-p38a interface. MK5 is 
in grey ribbon, while p38a is in yellow ribbon. Dashed lines indicate interactions between 
amino acids in MK5 and p38a. The percentages indicate the percentages of MD frames 
with the particular interacting atomic distance <3.5 A. (A) the activation segment of MK5; 
(B) the C-terminal end of MK5; (C) the P-loop of MK5; (D) the aH helix of MK5; 
(E) the regulatory phosphorylation region of MK5. 
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Calculations of the Electrostatic Potential Surfaces (EPS) indicated that the C-terminal segment 
with the p38 docking site and the regulatory phosphorylation region of MK5 are strongly 
electropositive and the other regions of MK5 facing p38a are mainly electropositive (Figure 11). 
However, the ATP binding pocket of MK5 located between the N-lobe and C-lobe is electronegative 
(Figure IIC). The regions of MK5 facing the solvent are partially electronegative and partially 
electropositive. p38a has a solvent-exposed side that is mainly electronegative, while the side facing 
MK5 is mixed electronegative and electropositive with highly electronegative areas in the regions 
interacting with MK5, the P-loop, helix aH, regulatory phosphorylation site, and p38 docking domain 
(Figure 11). This observation indicates that the charge distribution on the surface of MK5 and p38a is 
important for orienting the kinases relative to each other for a proper binding and that the binding is 
initiated by electrostatic forces. 

After optimization and equilibration the complex was simulated for 100 ns. The RMSF during MD 
(Figure 5) showed that the fluctuations of the flexible loops (including the P-loop) during MD were 
decreased compared with the simulations with MK5 alone and the inhibitor complexes. However, 
monitoring the atomic distance between Lys5 1 and Glu62 at the active site showed that this distance 
varied a great deal during the MD (Figure 7). A similar trend was also seen for the following atomic 
distances: Glu62-Phel70, Lys51-Ser33, Aspl69-Ser33 and Aspl69-ne32 (data not shown). These 
atomic distances were very stable during the MD with the free MK5 model. Therefore, it seems like 
p38a binding not only reduces the structural motions of flexible MK5 loops at the interaction surface, 
but also destabilizes intramolecular interactions between amino acids in the region of the active site, 
which may be important for MK5 activation and the binding of substrates to MK5. 

In general, the MK5-p38a complex was structurally very stable throughout the 100 ns MD. During 
the entire MD, amino acids in the P-loop, activation segment, aH helix, the regulatory phosphorylation 
region and the C-terminal were binding to amino acids of p38a (Figure 12). Table 1 show that amino 
acids in all these regions of MK5 formed stable interactions with p38a. In the activation segment the 
region Vall90-Glnl93 seemed important for binding and formed hydrogen bonds with amino acids in 
p38a. As previously suggested [14], the C-terminal segment was important for binding to p38a, and 
both Asn357 and Arg363 of the suggested p38 docking site interacted with p38a. During MD, Asn357 
formed a hydrogen bond with Aspl61, while Arg363 forms a salt bridge with Asp316 and a hydrogen 
bond with Tyrl32 of p38a. Outside the suggested p38 docking site in the C-terminal segment, Gln343 
formed a hydrogen bond with AsnllS of p38a that was observed in 99.4% of the sampled frames, 
while Lys346 had a salt bridge with Aspll2 that was observed in 63.8% of the frames. In addition 
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Lys350 had strong interactions with both Glul60 and Aspl61 (Table 1). In the P-loop, the backbone of 
Gly31 and Gly34 interacted with the backbone of Lysl5 and Glul2, respectively (Table 1). The 
backbone interaction between GlySl (P-loop) and Lysl5 in p38a was seen in 100% of the frames 
(Table 1). In aH helix, Asp322 formed salt bridges with Arg57, while Gln331, Gln335, Gln336 and 
Asn339 also formed stable interactions with p38a. In the regulatory phosphorylation region of MK5, 
Tyr242, Lys244 and Arg248 formed very stable interactions with p38a and were the most important 
amino acids for p38a interactions in this part of MK5. 

A combined molecular modelling and site-directed mutagenesis study suggested that Leul56 of 
p38a is important for directing the MK5-p38a complex to subcellular compartments [57]. When 
Leul56 was mutated to Val, the complex was mainly localized in the cytoplasm, and it was suggested 
that Leul56 interacts with the NLS. Our model did not show any direct interaction between Leul56 of 
p38a and MK5 (Table 1). However, our model showed that mutating Leul56 into Val will have an 
effect on the 3D structure of p38a, such that the observed interaction of Glul60/Aspl61 of p38a with 
the region around Lys350 of MK5 (Table 1) is disturbed. Lys350 is a part of the NES (Figure 13), and 
compared to the wild type p38a, the Leul56Val mutant may interact differently with the NES. 

Figure 13. A close up of the nuclear export signals (NES) and nuclear localization signals 
(NLS) of MK5. The MK5-p38a complex indicates that NES and NLS are masked by the 
binding of p38a. MK5 in grey ribbon and p38a golden color ribbon. 




3. Experimental Section 

The MD simulations were performed with the parallel MD program Desmond [58]. The following 
MD simulations were performed: 

1. 100 nsMD of MK5 alone. 

2. 100 ns MD of MK5 in complex with a pyrrolopyrimidone-based MK5 inhibitor (Chembel 
compound id: CHEMBL461 139). The inhibitor was in a neutral form. 
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3. 200 ns MD of MK5 in complex with the same inhibitor. The inhibitor was in a neutral form. 

4. 100 ns of MK5 in complex with p38a. 

3.1. Construction of Initial Complexes 

In a previous paper, we constructed homology models MK5 based on different templates using the 
ICM- software [45]. The models were tested by docking a library of known small molecular inhibitors 
and decoys. The test docking showed that the MK5 model based on the X-ray structure of MK3 in 
complex with an inhibitor (PDB id: 3FHR) [44] as template was capable of significantly splitting 
between binders and decoys. The MK3 X-ray structure is in an active kinase conformation and thereby 
also the MK5 model. The model consisted of the MK5 residues 9 to 348. The 3FHR-based model was 
used as a starting structure for MD simulations of MK5 alone and in complex with the inhibitor. The 
inhibitor was docked into the ATP binding pocket as described in our previous study [45]. In the 
previous study, we also constructed MK5 models based on the X-ray structure complex of MK2 and 
p38a (PDB id: 20ZA) [38]. In the present study, this model was used to obtain a starting complex of 
MK5 and p38 using the 20ZA complex as guide for the p38a docking. This model consisted of the 
MK5 residues 7 to 367. 

3.2. Calculation of Electrostatic Potentials 

The EPS of the MK5 and p38 were calculated with the Rapid Exact-Boundary Electrostatics 
(REBEL) method as implemented in the ICM program [59]. This method solves the Poisson equation 
for a molecule using a boundary element algorithm and generates a 3D surface skin model coloured by 
electrostatic potential. The potential scale values were set to from -5 to +5 kcal/e.u. charge units. At 
these scale values, the surface is coloured according the potentials: red potentials <-5, blue potentials 
>+5. A dielectric constant of 4 was used in the calculations. 

3.3. Solvation and Refinements 

The complexes were optimized with the protein preparation wizard in Maestro V9.1 [60] by 
assigning bond orders, adding hydrogen and correcting wrong bond types. The molecular systems were 
neutralized by adding chloride. In Desmond, the volume of space in which the simulation takes place 
(the global cell) was divided into regular 3D simulation boxes oflOAx lOAx 10 A. Each box was 
assigned to a single Desmond process. These boxes constituted the total simulation space with a total 
volume of 90 A X 90 A X 90 A. The molecular systems were solvated by SPC (simple point charge) 
Orthorhombic water box. 

All complexes were refined and energy optimized using the default quick relaxation protocol of the 
Desmond program (D. E. Shaw Research, New York, NY, USA). In short, two rounds of steepest 
descent minimization were performed with a maximum of 2000 steps with and without restraints 
(force constant of 50 kcal/mol/A on all solute atoms). These minimizations were followed by a series 
of four short MD simulations: (1) 12 ps MD simulation at a temperature of 10 K in the Berendsen 
NVT ensemble (constant number of particles, volume, and temperature) with solute heavy atoms 
restrained (force constant of 50 kcal/mol/A); (2) Twelve ps simulation was performed at 10 K with the 
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same restraints as in (1), but with the Berendsen NPT ensemble (constant number of particles, 
pressure, and temperature); (3) Using the same restraint as in A and B, a 12 ps simulation was 
performed in which the temperature was raised to 300 K using the Berendsen NPT ensemble; 
(4) Finally, a 24 ps simulation at 300 K using the Berendsen NPT ensemble without restraints 
was used. 

3.4. The MD Simulations 

After the refinements and equilibration steps, MD simulations were performed for all the molecular 
systems using the OPLS 2005 force field [61]. The pressure was kept constant at one bar and the 
temperature was kept at 300 K, using the Nose-Hoover chain [62] and Martyna-Tobias-Klein 
methods [63]. The short-range and long-range Columbic interactions were calculated with a cut off 
radius of 9 A and with the Smooth particle mesh method [64] (Ewald tolerance: 1.0 x 10"^). The 
M- SHAKE algorithm was used to constrain bonds containing hydrogen atoms. All the simulations 
used a multistep RESPA integrator [65] with 2.0 fs time step for bonded interaction and short range 
non-bonded interactions, and 6.0 fs time step for the long-range non-bonded interactions. During MD, 
the energies were recorded every 1.2 ps, while the coordinates were recorded every 4.8 ps. 

3.5. Analysing the Simulations 

The simulations were analysed using the Maestro program of the Schrodinger suite of programs. 
The Root Mean Square Deviations (RMSD) from the initial structure and the root mean square 
fiuctuations (RMSF) were calculated during all MD simulations. All atoms were included in the 
calculations of the RMSF. Atomic distances important for ligand binding and structural changes in 
MK5 were calculated for all the simulations. 

4. Conclusions 

In the present study, homology models were used to explore the molecular interactions of MK5 
using MD simulations. The homology models have previously been shown to select known MK5 
binders in front of decoys during docking [45]. Monitoring energies and RMSD during simulations 
indicated that energetically and structurally stable complexes were formed. These observations 
indicate that the homology models were robust and of quality to be used for structural predictions. 

The MDs with MK5 alone and with the inhibitor were started from a MK5 conformation in which 
Lys5 1 in aC and Glu62 in P3 were in close proximity, which is consistent with an activated protein 
kinase conformation. This conformation was very stable during 100 ns MD of MK5 alone. During MD 
simulations with the inhibitor, the atomic distance between Lys5 1 and Glu62 was also very stable, but 
the inhibitor binding disrupted the hydrogen bonding network between other amino acids at the 
binding site and thereby influenced their molecular motions. The inhibitor binding pose was stable 
both during the 100 ns MD and the 200 ns MD. The pyrrolopyrimidine ring system interacted with a 
polar and charged region of MK5, including Lys51, while the pyridine ring interacted within a 
hydrophobic region. In addition, water molecules were forming hydrogen-bonding interactions with 
the inhibitor during MD (Figure 8). The simulations with the inhibitor indicated that the inhibitor binds 
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to a MK5 conformation consistent with an active protein kinase and thereby functions as an inhibitor 
by occupying the ATP binding site. The inhibitor binding affects the hydrogen bonding network 
between active site amino acids, which then increases the atomic distance between Lys51 and Glu62 
compared with a fully active protein kinase structure. 

Calculations of the ESP of MK5 and p38a showed that the interaction surface of MK5 is strongly 
electropositive, while the interaction surface of p38a is mainly electronegative, indicating that 
electrostatic interactions are important for recognition and binding. The MK5-p38a complex shows 
that the binding of p38a is masking both the NLS and NES of MK5 and that may explain why 
MK5-p38a complexes are found exclusively in the nucleus [57], which is in agreement with previous 
suggestions that the NLS is overlapping with the p38 docking site of MK5 [14,57]. The present study 
shows that p38a can form strong interactions with MK5. During 100 ns MD, both stable hydrogen 
bonds and salt bridges were observed between the proteins. During the entire MD, amino acids in the 
P-loop, activation segment, aH helix, the regulatory phosphorylation region and the C-terminal were 
binding to amino acids of p38a. In the design of small molecular compounds that interfere with the 
binding of p38a to MK5, it may therefore also be important to target central amino acids outside the 
p38 docking site of MK5 (Table 1) as well as amino acids within the p38 docking site. The present 
study contributes with increased insight into the structure of MK5 and its interactions with molecular 
partners. The information can be used in the design of selective compounds interfering with MK5 
activity which would be useful in the elucidation of the exact biological role of MK5. 

Acknowledgments 

Inger Lindin is enrolled in BIOSTRUCT, the national PhD school in structural biology, and the 
support from BIOSTRUCT is acknowledged. The work was also supported by CPU-hours from 
NOTUR — The Norwegian Metacentre for Computational Science. 

Conflicts of Interest 

The authors declare no conflict of interest. 
References 

1. Cargnello, M.; Roux, P.P. Activation and function of the MAPKs and their substrates, the 
MAPK-activated protein kinases. Microbiol Mol Biol Rev. 2011, 75, 50-83. 

2. Moens, U.; Kostenko, S. Structure and function of MK5/PRAK: The loner among the 
mitogen-activated protein kinase-activated protein kinases. Biol Chem. 2013, 394, 1 1 15-1 132. 

3. Kostenko, S.; Shiryaev, A.; Dumitriu, G.; Gerits, N.; Moens, U. Cross-talk between protein kinase 
A and the MAPK-activated protein kinases RSKl and MK5. Recept. Signal Trans duct. Res. 
2011, 57, 1-9. 

4. Gaestel, M. MAPKAP kinases-MKs — Two's company, three's a crowd. Nat. Rev. Mol Cell Biol. 
2006, 7, 120-130. 

5. Kyriakis, J.M.; Avruch, J. Mammalian MAPK signal transduction pathways activated by stress 
and inflammation: A 10-year update. Physiol. Rev. 2012, 92, 689-737. 



Int. J. Mol Set 2014, 15 



4899 



6. Setemes, O.M.; Mikalsen, T.; Johansen, B.; Michaelsen, E.; Armstrong, C.G.; Morrice, N.A.; 
Turgeon, B.; Meloche, S.; Moens, U.; Keyse, S.M. Activation of MK5/PRAK by the atypical 
MAP kinase ERK3 defines a novel signal transduction pathway. EMBO J. 2004, 25, 4780^791 . 

7. Schumacher, S.; Laass, K.; Kant, S.; Shi, Y.; Visel, A.; Gruber, A.D.; Kotlyarov, A.; Gaestel, M. 
Scaffolding by ERK3 regulates MK5 in development. EMBO J. 2004, 23, 4770^779. 

8. Ronkina, N.; Kotlyarov, A.; Gaestel, M. MK2 and MK3 — ^A pair of isoenzymes? Front. Biosci. 
2008, 13, 5511-5521. 

9. Shiryaev, A.; Moens, U. Mitogen-activated protein kinase p38 and MK2, MK3 and MK5: 
Menage a trois or menage a quatre? Cell. Signal. 2010, 22, 1 185-1 192. 

10. Ni, H.; Wang, X.S.; Diener, K.; Yao, Z. MAPKAPK5, a novel mitogen-activated protein kinase 
(MAPK)-activated protein kinase, is a substrate of the extracellular-regulated kinase (ERK) and 
p38 kinase. Biochem. Biophys. Res. Commun. 1998, 243, 492^96. 

11. New, L.; Jiang, Y.; Zhao, M.; Liu, K.; Zhu, W.; Flood, L.J.; Kato, Y.; Parry, G.C.; Han, J. PRAK, 
a novel protein kinase regulated by the p38 MAP kinase. EMBO J. 1998, 1 7, 3372-3384. 

12. Gerits, N.; Mikalsen, T.; Kostenko, S.; Shiryaev, A.; Johannessen, M.; Moens, U. Modulation of 
F-actin rearrangement by the cyclic AMP/cAMP-dependent protein kinase (PKA) pathway is 
mediated by MAPK- activated protein kinase 5 and requires PKA-induced nuclear export of MK5. 
J. Biol Chem. 2007, 282, 37232-37243. 

13. Anderson, D.R.; Meyers, M.J.; Vernier, W.F.; Mahoney, M.W.; Kurumbail, R.G.; Caspers, N.; 
Poda, G.I.; Schindler, J.F.; Reitz, D.B.; Mourey, R.J. Pyrrolopyridine inhibitors of mitogen-activated 
protein kinase-activated protein kinase 2 (MK-2). J. Med. Chem. 2007, 50, 2647-2654. 

14. Setemes, O.M.; Johansen, B.; Hegge, B.; Johannessen, M.; Keyse, S.M.; Moens, U. Both binding 
and activation of p38 mitogen-activated protein kinase (MAPK) play essential roles in regulation 
of the nucleocytoplasmic distribution of MAPK-activated protein kinase 5 by cellular stress. 
Mol. Cell. Biol. 2002, 22, 6931-6945. 

15. New, L.; Jiang, Y.; Han, J. Regulation of PRAK subcellular location by p38 MAP kinases. 
Mol. Biol. Cell 2003, 14, 2603-2616. 

16. Kostenko, S.; Shiryaev, A.; Gerits, N.; Dumitriu, G.; Klenow, H.; Johannessen, M.; Moens, U. 
Serine residue 115 of MAPK-activated protein kinase MK5 is crucial for its PKA-regulated 
nuclear export and biological function. Cell. Mol. Life Sci. 2011, 68, 847-862. 

17. Kant, S.; Schumacher, S.; Singh, M.K.; Kispert, A.; Kotlyarov, A.; Gaestel, M. Characterization 
of the atypical MAPK ERK4 and its activation of the MAPK-activated protein kinase MK5. 
J. Biol Chem. 2006, 281, 35511-35519. 

18. Hansen, C.A.; Bartek, J.; Jensen, S. A functional link between the human cell cycle-regulatory 
phosphatase Cdcl4A and the atypical mitogen-activated kinase Erk3. Cell Cycle 2008, 7, 325-334. 

19. Deleris, P.; Rousseau, J.; Coulombe, P.; Rodier, G.; Tanguay, P.L.; Meloche, S. Activation loop 
phosphorylation of the atypical MAP kinases ERK3 and ERK4 is required for binding, activation 
and cytoplasmic relocalization of MK5. Cell Physiol 2008, 277, 778-788. 

20. Sohal, D.; Yeatts, A.; Ye, K.; Pellagatti, A.; Zhou, L.; Pahanish, P.; Mo, Y.; Bhagat, T.; 
Mariadason, J.; Boultwood, J.; et al Meta-analysis of microarray studies reveals a novel 
hematopoietic progenitor cell signature and demonstrates feasibility of inter-platform data 
integration. PLoS One 2008, 3, e2965. 



Int. J. Mol Set 2014, 15 



4900 



21. Perander, M.; Keyse, S.M.; Seternes, O.M. Does MK5 reconcile classical and atypical MAP 
kinases? Front BioscL 2008, 75, 4617^624. 

22. Gerits, N.; Shiryaev, A.; Kostenko, S.; Klenow, H.; Shiryaeva, O.; Johannessen, M.; Moens, U. 
The transcriptional regulation and cell-specific expression of the MAPK-activated protein kinase 
MK5. Cell Mol Biol Lett. 2009, 14, 548-574. 

23. Tak, H.; Jang, E.; Kim, S.B.; Park, J.; Suk, J.; Yoon, Y.S.; Ahn, J.K.; Lee, J.H.; Joe, CO. 
14-3-3epsilon inhibits MK5 -mediated cell migration by disrupting F-actin polymerization. 
Cell Signal 2007, 79, 2379-2387. 

24. Kostenko, S.; Johannessen, M.; Moens, U. PKA-induced F-actin rearrangement requires 
phosphorylation of Hsp27 by the MAPKAP kinase MK5. Cell Signal 2009, 27, 712-718. 

25. Kostenko, S.; Dumitriu, G.; Laegreid, K.J.; Moens, U. Physiological roles of mitogen-activated- 
protein-kinase-activated p38-regulated/activated protein kinase. World J. Biol. Chem. 2011, 2, 
73-89. 

26. Stohr, N.; Kohn, M.; Lederer, M.; Glass, M.; Reinke, C; Singer, R.H.; Huttelmaier, S. IGF2BP1 
promotes cell migration by regulating MK5 and PTEN signaling. Genes Dev. 2012, 26, 176-189. 

27. Yoshizuka, N.; Chen, R.M.; Xu, Z.; Liao, R.; Hong, L.; Hu, W.Y.; Yu, G.; Han, J.; Chen, L.; 
Sun, P. A novel function of p38-regulated/activated kinase in endothelial cell migration and tumor 
angiogenesis. Mol Cell Biol 2012, 52, 606-618. 

28. Yoshizuka, N.; Lai, M.; Liao, R.; Cook, R.; Xiao, C; Han, J.; Sun, P. PRAK suppresses 
oncogenic ras-induced hematopoietic cancer development by antagonizing the JNK pathway. 
Mol Cancer Res. 2012, 10, 810-820. 

29. Kress, T.R.; Cannell, LG.; Brenkman, A.B.; Samans, B.; Gaestel, M.; Roepman, P.; 
Burgering, B.M.; Bushell, M.; Rosenwald, A.; Filers, M. The MK5/PRAK kinase and Myc form a 
negative feedback loop that is disrupted during colorectal tumorigenesis. Mol Cell 2011, 41, 445^57. 

30. Tanguay, P.L.; Rodier, G.; Meloche, S. C-terminal domain phosphorylation of ERK3 controlled 
by Cdkl and Cdcl4 regulates its stability in mitosis. Biochem. J. 2010, 428, 103-1 11. 

31. Coulombe, P.; Rodier, G.; Pelletier, S.; Pellerin, J.; Meloche, S. Rapid turnover of extracellular 
signal-regulated kinase 3 by the ubiquitin-proteasome pathway defines a novel paradigm of 
mitogen-activated protein kinase regulation during cellular differentiation. Mol. Cell. Biol. 2003, 
23, 4542-4558. 

32. Chow, K.T.; Timblin, G.A.; McWhirter, S.M.; Schlissel, M.S. MK5 activates Rag transcription 
via Foxol in developing B cells. /. Exp. Med. 2013, 210, 1621-1634. 

33. Koh, G.C.; Schreiber, M.F.; Bautista, R.; Maude, R.R.; Dunachie, S.; Limmathurotsakul, D.; 
Day, N.P.; Dougan, G.; Peacock, S.J. Host responses to melioidosis and tuberculosis are both 
dominated by interferon-mediated signaling. PLoS One 2013, 8, e54961. 

34. Gerits, N.; van Belle, W.; Moens, U. Transgenic mice expressing constitutive active MAPKAPK5 
display gender-dependent differences in exploration and activity. Behav. Brain Funct. 2007, 3, 58. 

35. Brand, F.; Schumacher, S.; Kant, S.; Menon, M.B.; Simon, R.; Turgeon, B.; Britsch, S.; Meloche, S.; 
Gaestel, M.; Kotlyarov, A. The extracellular signal-regulated kinase 3 (mitogen-activated protein 
kinase 6 [MAPK6])-MAPK-activated protein kinase 5 signaling complex regulates septin 
function and dendrite morphology. Mol Cell Biol 2012, 32, lAei-lAl'i. 



Int. J. Mol Set 2014, 15 



4901 



36. Sun, P.; Yoshizuka, N.; New, L.; Moser, B.A.; Li, Y.; Liao, R.; Xie, C; Chen, J.; Deng, Q.; 
Yamout, M.; et aL PRAK is essential for ras-induced senescence and tumor suppression. Cell 
2007, 128, 295-308. 

37. Shi, Y.; Kotlyarov, A.; Laabeta, K.; Gruber, A.D.; Butt, E.; Marcus, K.; Meyer, H.E.; Friedrich, A.; 
Volk, H.D.; Gaestel, M. EUmination of protein kinase MK5/PRAK activity by targeted 
homologous recombination. Mol Cell. Biol. 2003, 23, 7732-7741. 

38. White, A.; Pargellis, C.A.; Studts, J.M.; Wemeburg, E.G.; Farmer, B.T., 2nd. Molecular basis of 
MAPK-activated protein kinase 2:p38 assembly. Proc. Natl. Acad. Sci. USA 2007, 104, 6353-6358. 

39. Underwood, K.W.; Parris, K.D.; Federico, E.; Mosyak, L.; Czerwinski, R.M.; Shane, T.; 
Taylor, M.; Svenson, K.; Liu, Y.; Hsiao, C.L.; et aL Catalytically active MAP KAP kinase 2 
structures in complex with staurosporine and ADP reveal differences with the autoinhibited 
enzyme. Structure 2003, 77, 627-636. 

40. Revesz, L.; Schlapbach, A.; Aichholz, R.; Feifel, R.; Hawtin, S.; Heng, R.; Hiestand, P.; Jahnke, W.; 
Koch, G.; Kroemer, M.; et al. In vivo and in vitro SAR of tetracyclic MAPKAP-K2 (MK2) 
inhibitors. Part I. Bioorg. Med. Chem. Lett. 2010, 20, 4715-4718. 

41. Revesz, L.; Schlapbach, A.; Aichholz, R.; Dawson, J.; Feifel, R.; Hawtin, S.; Littlewood-Evans, A.; 
Koch, G.; Kroemer, M.; Mobitz, H.; et al. In vivo and in vitro SAR of tetracyclic MAPKAP-K2 
(MK2) inhibitors. Part II. Bioorg. Med Chem. Lett. 2010, 20, 4719-4723. 

42. Fujino, A.; Fukushima, K.; Kubota, T.; Kosugi, T.; Takimoto-Kamimura, M. Crystal structure of 
human cyclin-dependent kinase-2 complex with MK2 inhibitor TEI-I01800: Insight into the 
selectivity. J. Synchrotron Radiat. 2013, 20 (Pt. 6), 905-909. 

43. Barf, T.; Kaptein, A.; de Wilde, S.; van der Heijden, R.; van Someren, R.; Demont, D.; 
Schultz-Fademrecht, C; Versteegh, J.; van Zeeland, M.; Seegers, N.; et al. Structure-based lead 
identification of ATP-competitive MK2 inhibitors. Bioorg. Med. Chem. Lett. 2011, 21, 3818-3822. 

44. Cheng, R.; Felicetti, B.; Palan, S.; Toogood-Johnson, I.; Scheich, C; Barker, J.; Whittaker, M.; 
Hesterkamp, T. High-resolution crystal structure of human Mapkap kinase 3 in complex with a 
high affinity ligand. Protein Sci. 2010, 19, 168-173. 

45. Lindin, I.; Wuxiuer, Y.; Kufareva, I.; Abagyan, R.; Moens, U.; Sylte, I.; Ravna, A.W. Homology 
modeling and ligand docking of Mitogen-activated protein kinase-activated protein kinase 5 
(MK5). Theor. Biol. Med Model 2013, 10, 56. 

46. Emrick, M.A.; Lee, T.; Starkey, P.J.; Mumby, M.C.; Resing, K.A.; Ahn, N.G. The gatekeeper 
residue controls autoactivation of ERK2 via a pathway of intramolecular connectivity. Proc. Natl. 
Acad Sci. USA 2006, 103, 18101-18106. 

47. Roskoski, R., Jr. ERKl/2 MAP kinases: Structure, function, and regulation. Pharmacol. Res. 
2012, 66, 105-143. 

48. Nicolau, N., Jr.; Giuliatti, S. Modeling and molecular dynamics of the intrinsically disordered e7 
proteins from high- and low-risk types of human papillomavirus. J. Mol. Model. 2013, 19, 
4025^037. 

49. Cashman, D.J.; Ortega, D.R.; Zhulin, LB.; Baudry, J. Homology modeling of the CheW coupling 
protein of the chemotaxis signaling complex. PLoS One 2013, 8, e70705. 

50. Wang, Z.; Li, S.; Sun, L.; Fan, J.; Liu, Z. Comparative analyses of lipoprotein lipase, hepatic lipase, 
and endothelial lipase, and their binding properties with known inhibitors. PLoS One 2013, 8, e72146. 



Int. J. Mol Set 2014, 15 



4902 



51. Gabrielsen, M.; Ravna, A.W.; Kristiansen, K.; Sylte, I. Substrate binding and translocation of the 
serotonin transporter studied by docking and molecular dynamics simulations. J. Mol. Model. 
2012, 18, 1073-1085. 

52. Rahman, K.S.; Cui, G.; Harvey, S.C.; McCarty, N.A. Modeling the conformational changes 
underlying channel opening in CFTR. PLoS One 2013, 5, e74574. 

53. Yuan, S.; Wu, R.; Latek, D.; Trzaskowski, B.; Filipek, S. Lipid receptor SlP(l) activation scheme 
concluded from microsecond all-atom molecular dynamics simulations. PLoS Comput. Biol. 
2013, P,el003261. 

54. Yuan, S.; Vogel, H.; Filipek, S. The role of water and sodium ions in the activation of the 
mu-opioid receptor, ^wg^w. Chem. Int. Ed. Engl. 2013, 52, 10112-10115. 

55. Arumugam, K.; Crouzy, S.; Chevigne, A.; Seguin-Devaux, C; Schmit, J.C. Structure prediction 
of GPCRs using piecewise homologs and application to the human CCR5 chemokine receptor: 
Validation through agonist and antagonist docking. J. Biomol. Struct. Dyn. 2013, 
doi:10.1080/07391 102.2013.817952. 

56. Schlapbach, A.; Feifel, R.; Hawtin, S.; Heng, R.; Koch, G.; Moebitz, H.; Revesz, L.; Scheufler, C; 
Velcicky, J.; Waelchli, R.; et al. Pyrrolo-pyrimidones: A novel class of MK2 inhibitors with 
potent cellular activity. Bioorg. Med. Chem. Lett. 2008, 18, 6142-6146. 

57. Li, Q.; Zhang, N.; Zhang, D.; Wang, Y.; Lin, T.; Zhou, H.; Ye, Z.; Zhang, F.; Lin, S.C.; Han, J. 
Determinants that control the distinct subcellular localization of p38alpha-PRAK and 
p38beta-PRAK complexes. J. Biol. Chem. 2008, 283, 11014-11023. 

58. Shan, Y.; Seeliger, M.A.; Eastwood, M.P.; Frank, F.; Xu, H.; Jensen, M.O.; Dror, R.O.; Kuriyan, J.; 
Shaw, D.E. A conserved protonation-dependent switch controls drug binding in the Abl kinase. 
Proc. Natl. Acad Sci. USA 2009, 106, 139-144. 

59. Totrov, M.; Abagyan, R. Rapid boundary element solvation electrostatics calculations in folding 
simulations: Successful folding of a 23-residue peptide. Biopolymers 2001, 60, 124-133. 

60. Maestro, 9.1; Schrodinger: LLC, New York, NY, USA, 2010. 

61. Jorgensen, J.H.; Swenson, J.M.; Tenover, F.C.; Barry, A.; Ferraro, M.J.; Murray, P.R.; Reller, L.B. 
Development of interpretive criteria and quality control limits for macrolide and clindamycin 
susceptibility testing of Streptococcus pneumoniae. J. Clin. Microbiol. 1996, 34, 2679-2684. 

62. Martyna, G.J.; Klein, M.L.; Tuckerman, M. Nose-Hoover chains: The canonical ensemble via 
continuous dynamics. /. Chem. Phys. 1992, 97, 2635. 

63. Martyna, G.J.; Tobias, D.J.; Klein, M.L. Constant pressure molecular dynamics algorithms. 
J. Chem. Phys. 1994, 101, Mil. 

64. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A smooth particle 
mesh Ewald method. Chem. Phys. 1995, 103, %511 . 

65. Humphreys, D.D.; Friesner, R.A.; Berne, B.J. A Multiple-time- step molecular dynamics 
algorithm for macromolecules. J. Chem. Phys. 1994, 98, 6885-6892. 



© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article 
distributed under the terms and conditions of the Creative Commons Attribution license 
(http ://creativecommons . org/licenses/by/3 . 0/) . 



